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We describe an experimental protocol for introducing spin-dependent lattice structure in a cold 
atomic fermi gas using lasers. It can be used to realize Hubbard models whose hopping parameters 
depend on spin and whose interaction strength can be controlled with an external magnetic field. We 
suggest that exotic superfluidities will arise in this framework. An especially interesting possibility is 
a class of states that support coexisting superfluid and normal components, even at zero temperature. 
The quantity of normal component varies with external parameters. We discuss some aspects of the 
quantum phase transition that arises at the point where it vanishes. 
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I. INTRODUCTION 



Cold atom systems can be used to explore impor- 
tant problems of condensed matter physics in new ways. 
For example, recent rapid development of the ultracold 
atomic gas in optical lattices 0, supported the obser- 
vation of the superfluid-Mott insulator transition in cold 
atomic bosons confined in an optical lattice 0, 0, IE E| In 
a wider context, methods of "engineering" various lattice 
model systems with bosonic and fermionic atoms have 
been proposed 0, H, 0, E3, EH E3, E3, E3, EH , 

open- 
ing prospects for exploring exotic new phases. In recent 
years achieving superfluidity in cold atomic fermions has 
become a major goa l, involving several experimental pro- 
grams [13 El M |U El Si IH S HI. Here we 
consider techniques that permit one to introduce both 
sorts of complexity in one system, with controlled band 
structures and interactions that depend on spin. They 
involve counterpropagating laser beams, that together 
generate a standing light wave which leads to different 
AC Stark shifts for the spin up and spin down com- 
ponents of alkali atoms (such as 40 K) in their ground 
state 0, Hi ■ We then discuss an interesting new phase 
of matter that might arise in this context, involving co- 
existence of normal Fermi liquid and superfluid compo- 
nents HHHm |H IS , and the quantum phase 
transition between this state and conventional Bardeen- 
Cooper-Schrieffer [33| (BCS) superfluidity. 



A. The setup of an optical lattice model 

Atoms in an off-resonant laser field exhibit a second 
order AC Stark shift of their ground state levels. This 
shift is proportional to the light intensity. The intensity 
may vary in space, for example, by forming a standing 
light wave from counterpropagating laser beams. For the 
center-of-mass motion of the atoms, spatially varying AC 
Stark shifts play the role of a conservative potential. In 
particular, a standing light wave leads to a periodic in- 
tensity pattern which results in a periodic potential, i.e. 
an optical lattice. By superimposing lattice beams from 
different directions and different intensity an effectively 
ID, 2D or 3D lattice can be built. 

Consider a cloud of cold fermionic atoms in 3D optical 
lattice EH- We assume that only two of the internal 
ground states participate in the dynamics, which we call 
spin up and down, | f) and | J,}. The corresponding 
effective Hamiltonian is 
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II. SPIN-DEPENDENT OPTICAL LATTICES 
FOR COLD FERMIONIC ATOMS 



In this section we describe how to realize, using 
cold atoms in an optical lattice, a rather general band 
Hubbard-type model, with separate, tunable effective 
masses and filling factors for the two spins 0, H E| ■ 



with i/v = i-^(x) fermionic field operators obeying the 
usual anticommutation relations. The various terms con- 
tributing to the Hamiltonian are illustrated in Fig. ^ 

The first line in (JJJ contains the the kinetic energy and 
the optical lattice potential V a (x) generated by the laser 
beams, obtained in second order perturbation theory for 
coupling of the ground state levels to the excited atomic 
states. In the simplest case of three orthogonal laser 
beams this potential has the form 
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FIG. 1: Spin-dependent lattice: two counterpropagating laser 
beams generate a standing light wave which leads to different 
AC Stark shifts for the spin up and spin down components of 
the ground state atoms (a =|, J.). As a result the tunneling 
matrix elements t a for the two spin components are different. 
The spin up and down components are coupled by a Raman 
laser with Rabi frequency Qq and detuning S. 



Here k = 2tt/X is the wave vector of the light, is 
proportional to the lattice beam intensity in the direc- 
tion t and the dynamic atomic polarizability at the laser 
frequency of the level (spin state) a. 

Note that in Eq. we have ignored cross AC Stark 
terms coupling two different spin components. These are 
negligible for the situation indicated in Fig. [2 where the 
AC Stark shifts are much smaller than the (bare) energy 
of the atomic levels, which can be achieved, e.g., by ap- 
plying a magnetic field to split the magnetic sublevels 
of the atoms. Furthermore, we have assumed in J2J that 
the AC Stark shifts of lattice beams in different directions 
can be added. Depending on the laser configuration there 
may be interference terms. However, for lattice beams in 
different directions with (slightly) shifted optical frequen- 
cies (obtained by modulating the lattice beams), these 
interference terms will average to zero. We have also ig- 
nored contributions from spontaneous emission. These 
can, for sufficiently large detunings of the lattice beams 
from the excited states, be made arbitrarily small. In- 
deed, for large detuning A ^> T from excited states with 
width r, the spontaneous emission rate scales as T/A 2 
while the lattice potential scales as 1/A, so spontaneous 
emission becomes, on the relevant dynamical timescale, 
negligible. 

A key element in JTJ is the assumption that the lattice 
potential V^(x) is spin dependent, as illustrated graph- 
ically in Fig. ^ Immediately below we will discuss in 
detail specific atomic and laser configurations that allow 
us, by varying laser parameters, to engineer this spin de- 
pendence. 

The second line in Eq. JTJ describes the coupling of 
the two spin states via a Raman transition with effective 
two-photon Rabi frequency firjj an d Raman detuning S. 
The Hamiltonian has been written in the frame, where 
the optical frequencies have been transformed away, so 



that only the detuning 6 appears in £[J. 

Finally, the last term in Q is the usual atom-atom in- 
teraction described by a pseudo potential with scattering 
length a s . In atomic experiments the sign and magnitude 
of the scattering length can be manipulated, for example, 
with an external magnetic field via Feshbach resonances 
pll I22I |23 | .The derivation of assumes explicitly that 
the scattering length is much smaller than the lattice 
spacing (given by A/2). In more general circumstances, 
the interaction will extend beyond nearest neighbors. 

For single atoms the energy eigenstates in the optical 
lattice are conveniently described by a band structure 
and Bloch wave functions. An appropriate superposi- 
tion of the Bloch wave functions forms a set of Wannier 
states which are localized at the various lattices sites. 
For typical experimental parameters the frequencies as- 
sociated with dynamics of cold atoms in the lattice are 
much smaller than the excitation energies in the optical 
potential (excitation to higher Bloch bands). Expanding 
the field operators in the set of Wannier functions of the 
lowest Bloch band, — J2i w <r( x ~ *-i) a ai, we obtain 

the Hubbard (single band) Hamiltonian 



+ "t H ( C W + h - c -) ~ u AAiW* ■ ( 3 ) 



In the Hamiltonian ©, h = —5/2 plays the role 
of "magnetic field" , t a is a spin dependent hopping 
term which follows from the spin dependent optical po- 
tential, and U is an onsite interaction. The depen- 
dence of the hopping amplitude t a and onsite interac- 
tion U on the optical lattice parameters is given by 
t a = E R (2/^F)(VojE R f/ 4 eM-^(Voa/E R ) 1 / 2 ) and 
U = E R ka a y/S/Tt(V 0<r I E R f/ 4 with E R = h 2 k 2 /2m the 
recoil frequency of the atoms, and k = 2*k /A the wave 
vector of the light. Thus the ratio of tunneling to onsite 
interaction can be controlled via the depth of the optical 
lattice. 



B. Tuning spin-dependent tunneling 

For alkali atoms, i.e. the atoms used in present cold 
fermi gas experiments, the difficulty in obtaining a spin 
dependent lattice arises from the s-wave character of the 
ground state. It implies that the AC Stark shift in- 
duced by far off-resonant driving fields is the same for 
all (hyperfine) ground state levels. Fortunately, however, 
heavy alkali atoms such as 40 K , have a larg e fine struc- 
ture splitting of the first excited state PJ|. This per- 
mits us to obtain a spin dependent optical potential by 
tuning the lattice lasers between the fine structure lev- 
els, still remaining far off resonance to suppress sponta- 
neous emission 8] . As illustrated in Fig. [2 the ground 
states n 2 5i/2 M = ±1/2 states couple to excited states 
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n 2 P-i/2,i/2 with right (respectively, left) circularly polar- 
ized light a according to the selection rules AM = ±1. 
This gives rise to a spin dependent AC Stark shift. 

For example, the AC Stark shift of the two ground 
states in circularly polarized light with amplitude £ is 



coefficients in the decomposition are sufficiently differ- 
ent. These spin dependent lattices were first proposed 
for quantum computing purposes by Jaksch et al. @ 
and recently implemented in experiments with cold Rb 
atoms by Bloch and collaborators [5| ■ 



_ 1^3/ 



|A»3/2^| : 



AE2 Sl/2M=+1/2 (<7 + ) 

A C 1 f~+\ Imi/2^1 2 i Ia»3/2£| 
Ab2 Si/2M= _ 1 / 2 (o- ) = 



fi(A-A fs ) 

, 1^3/2- 

hA fi(A-A fa ) ' 



(4) 
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keeping only the dominant contributions from the lowest 
lying excited states. Here the detuning from the 2 Pz/2 
state is denoted A, Af s is the fine structure splitting, and 
the dipole matrix elements for the transitions from the 
ground state are ^3/2,1/2- We see that for detunings be- 
tween the fine structure states < A < Af s the AC Stark 
shift AE2 Si/2M= _ 1 / 2 switches from positive to negative 
values, while AE2 Si , 2 m=+i/2 is always negative, giving 
rise to a strongly spin dependent lattice potential. We 
can vary the strength of the Stark shift by varying laser 
power. 

Thus, the laser configuration consisting of unbalanced 
right and left circularly polarized polarization compo- 
nents gives rise to a spin dependent lattice, where the 
resulting AC Stark shifts are a sum of the shifts in the 
first and second line of (0J, weighted according to their 
fraction of a + and o~ component. For linearly polarized 
light the AC Stark shifts of the two states are equal. 
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FIG. 2: a) Atomic level scheme of Alkali atoms. A a + polar- 
ized lattice beam is tuned between the two excited state fine 
structure components resulting in different AC Stark shifts for 
the ground state levels as explained in the context of Eq. @ • 
b) Two possible "spin up" and "spin down" states are illus- 
trated the case of 40 K atoms. 

It is easy to include hyperfine splitting of the atomic 
ground states Q . Consider an atomic ground state elec- 
tron which is coupled to a nuclear spin I. The resulting 
total angular momentum is F = \I ± 1/2 1, and the hy- 
perfine ground states have the form \F, Mp) = a\I, Mp + 
1/2) | 2 S 1/2 M = -l/2)+b\I,M F -l/2)\ 2 S 1/2 M = +1/2) 
where a and b are Clebsch-Gordan coefficients. Assuming 
that the Zeeman hyperfine states are split by a constant 
magnetic field (compare the discussion following Q), the 
AC Stark shift of this state will be the weighted sum 
AEpm f = \a\ 2 AE2 Si/2M= „ 1/2 + \b\ 2 AE2 Si/2M=+1/2 . 
Thus a pair of hyperfine ground state levels will have a 
spin dependent lattice provided that the Clebsch-Gordan 
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FIG. 3: AC Starkshift (in arbitrary units) of the 40 K hyperfine 
Mf = 9/2 (solid line) and Mf =7/2 (dashed line) states 
(compare Fig. [^Jd) as a function of the detuning A. The value 
of A = corresponds to the 2 P!/ 2 excited fine structure state, 
while A = Af s is the 2 P3/2 resonance. We note the strong 
spin dependence (i.e. dependence on the internal state) for 
detuning between the two fine structure states. 



This shows how to create a ID spin dependent lat- 
tice in counterpropagating laser beams of unbalanced 
polarization. The set-up is readily generalized to higher 
dimensions. For example, a 3D lattice is obtained by first 
applying a magnetic field to provide a Zeeman splitting 
of the hyperfine ground states (to suppress the cross AC 
Stark terms), and then applying three standing wave cr + 
polarized beams which are tilted by 45° relative to the 
magnetic field. Figs. [3] and ^summarize the correspond- 
ing results for the Mp = 9/2 and 7/2 hyperfine structure 
ground states (compare the level scheme in Fig. |2t>) ■ For 
these states the squares of the relevant Clebsch-Gordan 
coefficients (as defined above) are given by \a\ 2 = 1/9 
and |&| 2 = 8/9. Fig. plots the AC Stark shift of the 
M F = 9/2 (solid line) and M F = 7/2 states (dashed 
line) as a function of the laser detuning A. The detuning 
interval covers the region of the 2 Pi/2 and 2 P3/2 excited 
fine structure states, which are separated by Af s . We 
see that for a detuning A « 0.248 Af s the AC Starkshift 
of the Mp = 7/2 state has a zero while the shift of the 
Mp = 7/2 varies comparatively slowly. For detunings 
in the range 0.248 < A/A fs < 1 the AC Stark shift of 
both states Mp — 9/2,7/2 is negative, i.e. the minima 
of the optical wells V a (x) coincide but the wells have dif- 
ferent depths giving rise to a spin-dependent tunneling. 
The AC Starkshifts in this region between the two fine 
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structure states are plotted in Fig. 0Jl where again the 
Mp = 9/2 level is represented by the solid line, and the 
Mp = 7/2 level of 40 K is given by the dashed line. In 
Fig. 0Jd we give the ratio of the corresponding hopping 
matrix elements as a function of A in the same interval. 
Here we have chosen a light intensity so that we keep the 

,(0) 



depth of the optical potential K 



T=|JW>=9/2)/ 



V (for a 

given direction £ — 1,2,3) at a fixed given value Vb for 
the whole range of detunings, while ^ / |='|m f =7/2) i var i es - 
Thus the ratio of t\tjt\i varies. From Fig. 0Jd we see 
that the ratio of the tunneling elements can be varied by 
more than an order of magnitude when we approach the 
interference minimum at A 0.248 Af s . 

Besides the alkali atoms discussed above, experimental 
progress might soon allow the realization of quantum de- 
generate fermi gase s with alkaline earth atoms, e.g., with 
Sr atoms [35l l36(. Alkaline earth atoms have besides 
their singlet ground states longlived electronic excited 
triplet states. Applying an off-resonant laser field, these 
states will in general have quite different AC Stark shifts 
[37l | . Identifying the ground and metastable excited state 
with the spin-up and spin-down states, we thus can also 
have a natural realization of a Hubbard model with a 
spin-dependent interaction. 



C. Remarks on spin superposition by Rabi 
coupling 

The model we have been able (conceptually) to engi- 
neer contains a Rabi coupling that is not present in the 
conventional Hubbard model J3J). A few remarks about 
this issue are in order. 

Qq = corresponds to the case that the particle num- 
ber in each spin species is conserved separately. For 
£!o 7^ only total particle number conserved, while the 
relative particle number can vary. This point will be im- 
portant for our discussion of exotic phases. 

To treat the effect of the Rabi term theoretically, we 
should first diagonalize the quadratic part of the Hamil- 
tonian. This is best done in momentum space. In prin- 
ciple, we can always find a unitary transformation of the 
original fermion fields, 



c ff (k) =^(k) c CT (k), 
that transforms the Hamiltonian into 

H = ^?ak2^ k 5ak 
ker 



(5) 



(6) 



^ ' Kr 1 ..CT 4 ;ki..k 4 C^ lkl cj r2k2 C ( j 3 k3C CT4 k4 

with V ai .. a4 . M .. ki = US d (k 1 +k 2 -k 3 -k 4 )^ lT (k 1 ) x 
^ T2 ^(k 2 )'^' cr ^(k3)'^'J 4l .(k4) in d-spacial dimensions. The 
explicit form of the unitary transformation matrix (k) 
for each mode k can be found easily after some elemen- 
tary algebra, but is not essential to our discussion here. 
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FIG. 4: a) AC Starkshift (in arbitrary units) of the Mp = 9/2 
(solid line) and the Mp = 7/2 state (dashed line) in the de- 
tuning region right of the interference zero 0.248 < A/Af s 
of the Mp = 7/2 state, b) Ratio of the hopping matrix ele- 
ments for the t^=\u F =9/2),i/ti=\M P =7/2),e for a given spatial 
direction £ = 1,2,3 as a function of detuning A for a fixed 
VjS Mp=g t 2 K e = Vb with Vo = 5 (curve a), 10 (curve b), and 
20Er (curve c) in units of the recoil energy Er = h 2 k 2 /2m. 
The hopping matrix elements were obtained from a band 
structure calculation for the given depth of the optical po- 
tential. 



The important point is that the new, effective energy 
bands e a k remain spin dependent. The free parameters 
h, f2o, and t^s allow great freedom of tuning the new 
band mass ratio and the Fermi surface difference of the 
effective up and down fermions. Instead of working with 
the Hamiltonian © of the original fermions, we could 
just as well start with the Hamiltonian 10 of the effective 
fermions. The new interaction term will contain not only 
terms s-wave spin-singlet terms, but also other higher an- 
gular and spin terms. As far as our main interest, s-wave 
spin-singlet pairing, is concerned, the theoretical treat- 
ment does not differ substantially from the model Q 
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below, for which the Rabi term is transformed away. 



III. BP TO BCS TRANSITION 

Recently there has been considerable discussion of 
the possible existence of homogeneous zero-temperature 
phases wherein a superfluid condensation coexists with 
one or more Fermi surfaces, where the gap vanishes. 
We have in mind that these are full-fledged codimension 
1 Fermi surfaces (e.g., two-dimensional surfaces, which 
bound three-dimensional regions, in a three-dimensional 
system) where the gap vanishes. The roots of this idea 
go back to early work of Sarma (3^ ; the issue has arisen 
again, under several different names, in several new con- 
texts, including high-density QCD and (as emphasized 
here) cold atom systems. There are delicate issues of 
stability involved, which have been mishandled in much 
of the literature. We believe that a careful and cor- 
rect discussion is supplied in Ref. [3^|. We shall not 
repeat the analysis given there, but we rely on the conclu- 
sion: there are a number of physically interesting circum- 
stances in which two-component interacting Fermi sys- 
tems, of the general type discussed in the preceding sec- 
tion, can support a continuous "breached pairing" (BP) 
to BCS quantum phase transition (i.e., a "superfluid + 
normal— >superfluid" transition at zero temperature). 

Stability of the BP phase, which has coexisting, ho- 
mogeneous superfluid and normal components at zero- 
temperature, appears to require momentum-dependent 
interactions, same-species repulsion, or some other com- 
plication not present in the simplest Hubbard Hamilto- 
nian with two states defining a quasi-spin degree of free- 
dom. For the sake of simplicity and concreteness, how- 
ever, we shall work with this Hamiltonian. The qualita- 
tive, universal features of the phase transition ought not 
to depend on this idealization. Then, 

H = - t a (cltCvj + h.c.J +h^2 ( c |i c Ti " c h c li) 

- u J2 c \Ai c i- ic ^- ( 7 ) 



where as usual, c„i and c ai are fermion annihilation and 
creation operators at site i, with a =f, J, indicating two 
internal quantum states. 

The Hamiltonian Q appears to be a simplified ver- 
sion of ©, by assuming that the Rabi coupling term is 
transformed away through JSJ. An important difference 
is that the spin up and down internal states are coherent 
now. The two quasi-spin states could well be two hyper- 
fine spin levels for the case of cold atoms, as we discussed 
in detail above. 



A. Nature of quantum phase transition between 
BCS and BP states 

The general nature of the transition can be gleaned 
from Fig. [SJ The BP to BCS transition is continuous, 
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FIG. 5: Quantum phase transition seen from change of quasi- 
particle spectra. The critical value h = h c is determined by 
requiring that E% branch be of just one solution at zero en- 
ergy. For h > h c , one should imagine a three dimensional 
momentum space in which two surfaces, defined by Ok = 0~ 
and Ok = + , are gapless for the E2 branch, a) BCS super- 
fluid; b) Quantum critical state; c) "superfluid+normal" BP 
state. For convenience, t\ > £4, is always assumed throughout 
the paper. 

at least in mean field theory. To characterize its critical 
dynamics, we must use paramaeters that characterize he 
singular behavior — in particular, the appearance or dis- 
appearance of low-energy degrees of freedom. The magni- 
tude of the superconducting order parameter is evidently 
not a suitable order parameter, since it evolves smoothly 
(at the level of mean field theory) and is non- vanishing at 
both sides of the transition. Furthermore, there is no ev- 
ident change in symmetry at the transition. We propose 
that a good order parameter for the BP to BCS transition 
is the quasi-spin polarization (see Section [ill Dj l. 



B. Quasiparticle excitations 

We shall assume a uniform pairing, with order param- 
eter, A = U (el, cL) = const. To find the spectrum of low 
energy quasiparticle excitations, we diagonalize the mean 
field Hamiltonian 

^-2^ w kl a* -(2t,e k -/o + /i J 

k 

(8) 

where 
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and 



1=1. .d 



cos(fc;a) . 



(Here d = 3 for three dimensional space.) There are two 
physical branches of (Bogoliubov) quasiparticle excita- 
tions that have distinct property in a BP state. They 
are 



£ 1)2 (k) = (h + *_e k ) ± V(*+©k - /i) 2 + |A| 2 (9) 

where we have adopted the short-hand notation t± = 



ti±ti- We will assume tt > t \ and h > h 



)/'• 



By definition, /i x is the value of the magnetic field for 
which the spin up and down bands have the same Fermi 
surface. 

Fig.[5]shows how the two branches of excitations evolve 
with h. A quantum phase transition takes place at the 
critical value of h = h r with l40l 



-fit- + 2\A\y/t^tl 



(10) 



An important feature of the BP state is that the E 2 
branch of excitations crosses the zero energy axis at two 
two-dimensional surfaces in three-dimensional momen- 
tum space, defined to be Ok = 9~ and 9k = 9 + . Their 
values are 



9= 



(11) 

At the critical point h — h c , the two gapless "Fermi" 
surfaces merge into a single one, 



9" 



9^ 



9 C 



fit+ + h c t- 



(12) 



Excitation spectra similar to Fig. |3Jc) are also found 
in ferromagnetic metals. It was argued in |4lT | that su- 
perconductivity could coexist with fcrromagnetism. 



C. As a Lifshitz topological transition 

Lifshitz (topological) transitions ^2 can take place in 
metals and alloys at the motion of von Hove singularities 
of the electron density of states across Fermi surfaces (for 
review, see, e.g., Blanter et al. |43jp. The transition from 
BCS to BP actually falls into this universality class. 

To illustrate the topological nature of the transition, 
we start with the many-body wavefunctions for the BP 
and BCS states, 



\Vbcs) = U(uk + v k ct^_ H )\0), 

k 

i*bp> = n («w + «k4t c -ki) n 



(13) 



k:G k <e- k : e k e[e-,e+] 

x n K + ^Ck T cL ki )|0). (14) 
k : e k >e+ 



Here, «k and i?k are complex numbers satisfying |uk| 2 + 
|i>k| 2 = 1; their amplitudes are determined by diagonal- 
izing the mean field Hamiltonian (JSJ, 



l«k| 
\vu\ 




i+9 k - /i 



v/[i+0k-M] 2 + |A| 



(15) 



There is a manifold of degenerate states featuring an 
overall relative phase between the «k and «k factors, cor- 
responding to the broken U(l) c hargc- (Note that by our 
convention of h > h x , there are more particles in the 
spin I species than the f species. 

The occupation numbers of both spin up and down 
fcrmions are readily determined. For the BCS state, 



"-kT = "kl = \vk\' 

while for the BP state: 

if 9 k G [9-, 9+]: 
otherwise: 



(for all k) 



n k | = , n kl = 1 
n-kt = "ki = M 2 • 



(16) 



(17) 



In the BCS state the occupation numbers are equal, while 
for the BP state they differ. Since the BCS state is main- 
tained for a finite range of parameters, its continuous 
evolution to the BP state is non- analytic. This signals a 
(zero-temperature, quantum) phase transition. 

The phase transition associated with a change in topol- 
ogy of the Fermi sea. This kind of transition is gener- 
ally known as Lifshitz transition. As shown in Fig. El 
the quantum phase transition connects states of a simply 
connected Fermi sea and of isolated two regions. 




FIG. 6: An illustration of changing topology of the spin up 
Fermi sea. Shown are cross-sections at k z = n /2a in a 3D first 
Brillouin zone for (a) BCS and (b) BP; momenta units for 
k x , k y are n/a where a is the cubic lattice constant. As usual, 
the occupied states spread out over the entire Brillouin zone 
due to pairing interaction. The shading in the graph indicates 
that the region has a non-vanishing occupation number. The 
dotted lines indicate the orginal free Fermi surfaces of spin ups 
(inner curve) and spin downs (outer curve) . The solid lines in 
(b) are energy surfaces defined by Ok = Q~ and 6k = + ; 
For small A — > 0, they merge with the free Fermi surfaces 
(dotted lines), respectively. The phase transition occurs at 
the point in which a simply connected Fermi sea is isolated 
into two regions. 
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A transition between superfluid phases of different 
topology in momentum space is also known to occur in, 
for instance, 3 He-B phase. There, the superfluid velocity 
with respect to a container drives a quantum phase tran- 
sition from a fully gapped state to a state of (gapless) 
fermi surface I44ll4fil. 



D. The order parameter 

We now come to discuss the order parameter that char- 
acterizes the quantum phase transition. It is spin polar- 
ization, M = (5 z (x)) = n-f — n±. However, unlike a usual 
Landau-Ginzburg type theory, no obvious spontaneously 
broken symmetry is involved at the transition. The spin 
polarization density is given by 



M 



dQ N(Q) 



(18) 



where 



AT(e) = ^5(e-e k ) 



plays the role of (dimensionless) density of states. Near 
the critical point h = h c , both 9 _ and + approach 
C . So the spin polarization may be approximated to 
the lowest order in h — h c by 



M 



n(o c )(q- - e+) 



, , \/(ut- + ht + ) 2 - 4A 2 i U , x 

= -jv(e c )^- ^4 — , (19) 



2Ut 



-N(e c ) 



2Ut 



TM 



■y/(h-h c ){h + h c -2hx). 



The spin polarization may be thought being derived 
from the grand thermodynamical potential Q by viewing 
it as a functional of h: 



M = 



sn[h] 

Sh 



This formula can be inverted. Define Km as the mag- 
netic field for which the above equation has a prescribed 
value M. Then the quantum effective potential T[M] is 
defined (as a functional of M, not h) by the Legendre 
transformation 



T[M] = -KM + n[h] , 



with 



h = h 



fixed. The form of T[M] is determined by requiring that 
the equation of state, 



ar[M] 

dM 



(20) 



produce what is equivalent to Eq. I|19f) upon inverting the 
latter. This leads to 



r[M] 



t- , r AJhh 
— uM H — — 



+MAarcsinh 




\M A J 



(21) 



where Ma is a constant, Ma = AiV(9 c )/ v /FpfJ. In de- 
riving Eq. J2J, we have treated A as a magnetic field- 
independent parameter to effectively represent the cou- 
pling strength. The effective potential is not analytic in 
M; the appearance of \M\ is a consequence of selecting 
the physically stable solution of it. The physical origin 
is due to the presence of two gapless "Fermi" surfaces. 
For small M (and also M < in our particular case of 
ij > ti and h > h x ), we could expand T[M] in powers of 
M. Then, in the presence of external magnetic field h, 
the final effective potential reads 



V(M) = hM + T[M] 
= (h-h c )M 




(22) 



For h > h Cl V(M) always has a minimum at a nonzero 
M. From either Eq. (JTUJl or one can immediately 

verify the following power-law (scaling) relation for the 
spin polarization at mean field level (for h near h c ) 



M = 



N(e c )(t + A)? {h _ K)1/ ^ - lih>hc 



0. 



(*t*i) 3 



otherwise . 



(23) 



The results are sketched in Fig. [7] The behavior of 
spin polarization near the point h = h c suggests that the 
quantum phase transition is of second order. 




FIG. 7: Quantum phase transition shown by spin polarization 
M. For the BCS (BP) state, M = (M + 0). 
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E. Scaling theory of z — 2 gapless fermions 

At the critical point, the ^(k) branch of quasiparticlc 
excitations is gapless at a single 'Fermi surface' defined by 
0(k) = C . What makes the critical point special is that 
the spectrum disperses quadratically near the surface: 

-£ 2 (k) = -^, k = K + l (24) 

2771-2 

where K is a vector on the critical Fermi surface and 
1 a vector of (small) fluctuating momentum orthogonal 



to the Fermi surface. m* 2 is effective band mass that 
is implicitly determined from Eq. lO- The dispersion 
relation implies that the dynamical exponent is z = 2 
(at least at tree-graph level) in contrast with z = 1 in a 
nominal (metallic) Fermi liquid |46j| . 

The critical properties of the -E 2 (k) fermions appear to 
define a new universality class. To begin understanding 
it, let us consider the relevant operators. Consider a 
theory of gapless fermions described by an effective free 
action plus possible quartic interactions, 



1+ = J dtd^Kdl ^t(k)^ ft __L_ 



V-(k) 



gfdt]] [d^Kjdlj] ^ t (k 1 )^(k 2 )^(k 3 )V(k 4 ) ( 5 d (k 1 +k 2 -k3-k4) 

A 1/1 



j=l.A 

-j-lplplplfj + ij}' ijjljjlp . . . 



where tp is the gapless fermion field. 



(25) 



r 



One can perform a simple renormalization group analysis 
to find out how the interactions should scale. When we 
scale down the fluctuating momentum 1 by a factor s < 1, 

l-> si, 

the quasiparticle energy scales down with a dynamical 
exponent z = 2 as 

oj — > s 2 uj 

for a quadratic dispersion relationship. So, in geneneral, 
the momentum and time in the above action should scale 
as follows: 



dt 
d t 



„2< 



s z dt, d a ~ L K 



s 2 d t , 



l d -h 

l-al, 



K, dl 

with s — ► 0. 



s 1 dl, 



Note that the momentum K, being attached onto the 
critical Fermi surface, does not scale. Accordingly, the 
fermion field scales as 

The scaling dimension of potential interaction opera- 
tors can now be derived by straightforward power count- 
ing. We find that a generic four-fermion scattering op- 
erator is marginal. That is, 



9 



s°g 



(generic scatterings) 
s~ 1 g , (BCS or forward scatterings) 



This result contrasts dramatically with what is found 
in the standard renormalization group theory of a con- 
ventional Fermi liquid (see, for instance, the review by 



Shankar |46j1. In that context, only special momentum 
configurations are marginal. 

We next turn to the scaling property of the order pa- 
rameter correlation function. We now attempt to de- 
rive the critical theory of quantum fluctuations of the 
order parameter field M. While the static part of it is 
already derived in Eq. (|22|l . its dynamical part can be 
obtained by calculating the spin-spin correlation, which 
is just the fermion polarization bubble diagram, i.e., 
(^(k, iu>)S z (— k, — iu>)) — II(k, ioj) (in imaginary time 
formalism), where S z (k,iuj) is defined as the Fourier 
tansform of 



S z 



— ( ct ( 
2 V 



Ti 



at site i and time t. The quadratic part of the action for 
M is given by II(k, iui). Near the critical point, we find 
the quantum effective action of M takes the form, 



S[M] = d d xdT\^(d T M) 2 + ^-(WM) 2 

J L Z Z 

+ 2~ m2+ ! |m|3+ " 



(26) 



where 



J (^)* T (»- W », ,.J|L, (27, 



«(ej(f+A)i 



and 



Sn(k,icj) 
d(u 2 ) 



N 2 (e c )t+A 

an(k,iw) 
a(k 2 ) ' 



(28) 
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with k, lu — > (long wavelength limit). In deriving the ef- 
fective action, we have shifted the position of the minima 
of V(M) to eliminate the linear term of Eq. J22Jl- 

One recognizes that the effective theory of spin 
polarization field M bears superficial resemblance to 
Fisher's 0] itfi 3 famous field theory of the Yang-Lee edge 
singularity in a ferromagnetic Ising model. However, the 
two theories differ fundamentally in that the cubic inter- 
action term here, |M| 3 , is both non-analytic and real (as 
opposed to analytic and purely imaginary.) 

Further investigation of the new critical theory, includ- 
ing interactions between the low-energy fermion modes 
and the order parameter, is reserved for future work. 
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